Microwave Photon Counter Based on Josephson Junctions 
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We describe a microwave photon counter based on the current-biased Josephson junction. The 
junction is tuned to absorb single microwave photons from the incident field, after which it tunnels 
into a classically observable voltage state. Using two such detectors, we have performed a microwave 
version of the Hanbury Brown and Twiss experiment at 4 GHz and demonstrated a clear signature 
of photon bunching for a thermal source. The design is readily scalable to tens of parallelized 
junctions, a configuration that would allow number-resolved counting of microwave photons. 



PACS numbers: 42.50.-p,85.25.Pb,85.60.Gz 

The strong interaction of superconducting integrated 
circuits with microwave photons forms the basis of cir- 
cuit quantum electrodynamics (cQED), an attractive 
paradigm for scalable quantum computing and a test 
bed for quantum optics in the strong coupling regime 
[IHS]. While single-photon detectors are an important 
component of the standard quantum optics toolbox [7], 
they have not played a role in cQED, as there are no 
microwave photon counters currently available. In the 
search for a microwave component that will emulate the 
nonlinearity of traditional optical photon counters, it is 
natural to consider the Josephson junction - a nonlinear, 
nondissipative superconducting circuit element with tun- 
able microwave-frequency transitions. There has been 
recent interest in the development of on-chip detectors 
to probe the microwave radiation emitted by mesoscopic 
quantum circuits 8-14., including proposals to utilize 
Josephson circuits as photon counters in the context of 
cQED [51 H] . This study represents the first experimen- 
tal measurement of the coincidence counting statistics 
of microwave photons. We exploit microwave-induced 
transitions between discrete energy levels in a Josephson 
junction to realize an efficient, versatile microwave detec- 
tor. Moreover, we use a parallelized detector circuit to 
probe the counting statistics of both coherent and ther- 
mal microwave sources. 

A circuit diagram of the single-channel detector is 
shown in Fig. [TJa). The Josephson junction is current 
biased near the critical current Iq. The junction poten- 
tial energy is given by U(S) = — cos<5 — 7b<5, where 
/b is the bias current, i5 is the superconducting phase dif- 
ference across the junction, and $o = h/2e is the mag- 
netic flux quantum. For bias currents approaching the 
critical current, the local minima of the potential (char- 
acterized by barrier height AU and plasma frequency w p ) 
accommodate a handful of discrete energy levels. Reso- 
nant microwaves induce coherent oscillations between the 
ground state |0) and the excited state |1) of the junction 
(Fig. [ljb)). Moreover, each of these states can tunnel to 
the continuum at a rate r 0j i. Due to the exponential de- 
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FIG. 1. Microwave photon detector based on a Josephson 
junction, (a) Schematic diagram of the circuit. The junction 
is biased with a dc current, and microwaves are coupled to 
the junction via an on-chip capacitor. The voltage across the 
junction is read out by a room temperature preamplifier and 
comparator, (b) Junction potential energy landscape. The 
junction is initialized in the |0) state. An incident photon 
induces a transition to the |1) state, which rapidly tunnels 
to the continuum, (c) Calculated detection efficiency r\ as 
a function of applied microwave frequency and current bias. 
The simulation parameters are Iq = 210 fiA, C — 45 pF, Ti 
= 3 ns, 6JH./27T = 200 MHz, and tdet ~ 5 ns, chosen to match 
the parameters of our measured device, (d) Experimentally 
measured detection efficiency. The efficiency shows a maxi- 
mum of 0.7. The detection bandwidth is of order the Rabi 
frequency associated with the microwave drive. 



pendence of tunneling on barrier height, Tj is 2-3 orders 
of magnitude larger than Tq. 

Our experimental protocol involves pulsing the junc- 
tion bias for a finite interval (of order 10 ns) during which 
the transition frequency w 10 is near resonant with the in- 
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FIG. 2. (a) Experimental setup of microwave-frequency ver- 
sion of the HBT experiment. Coherent microwaves (upper left 
panel) are derived from a microwave generator whose output 
is coupled to the detector via a transmission line, with suitable 
attenuation. The thermal source (lower left panel) is realized 
by amplifying, filtering, and upconverting the Johnson noise 
of a room temperature resistor. A two-junction circuit (right 
panel) is used to probe the statistical nature of the microwave 
photon source; each junction has its own bias and readout, 
(b) Conditional switching probabilities Poo, Pio, -Poi, and P\\ 
versus relative measurement delay r for coherent microwaves 
(upper panel) and thermal microwaves (lower panel). Inset 
shows the timing diagram of the measurement pulses applied 
to the two junctions. 

cident photons. Absorption of a photon promotes the 
junction to the |1) state, which tunnels rapidly to the 
continuum, producing a large voltage pulse and a mea- 
surable "click" . We have calculated the time evolution of 
the junction density matrix in the presence of a classical 
coherent drive. We take into account T\ relaxation and 
tunneling to the continuum using standard open quan- 
tum systems tools [15]. To quantify the performance 
of the detector, we introduce the efficiency parameter 
V = -Pbri g ht(l - -Pdark), where Plight is the probabil- 
ity of a microwave-induced transition to the continuum 



during the measurement interval and Pdark is the prob- 
ability of a spurious dark count due to quantum tunnel- 
ing. We find that r\ is maximized for a bias point corre- 
sponding to AU/hiOp ~ 2 and for a measurement interval 
that is roughly equal to the Rabi period of the coherent 
drive; the detection bandwidth is of the order of the Rabi 
frequency u)r/2w. We note that for comparable to 
vacuum Rabi frequencies achieved in cQED experiments 
(wr/27t ~ 100 MHz), detection efficiencies in excess of 
0.7 are attained with junctions having extremely modest 
coherence times of order several ns. 

It is instructive to compare our Josephson photon 
counter to other Josephson junction-based detectors. Su- 
perconducting tunnel junction detectors are biased in the 
voltage state and incident photons (typically at near in- 
frared frequencies and beyond) break Cooper pairs to 
induce an enhanced subgap current [16J : in contrast, our 
junctions are biased in the supercurrent state, and inci- 
dent microwave photons induce a transition to the volt- 
age state. In other recent work, the bifurcation of a 
driven Josephson junction has been exploited to realize 
sensitive threshold detectors [T7] or amplifiers for field 
quadratures [18] . In contrast, our detector is strictly sen- 
sitive to the photon number: it registers a click whenever 
a nonzero number of photons is present, irrespective of 
that number. 

We have fabricated a double photon counter sample in 
which two Josephson junctions (labeled A and B) with 
independent biasing and readout are driven by a com- 
mon microwave source. The junctions were made in a 
standard A1-A10 X -A1 technology, with areas of 1000 zxm 2 . 
The detectors are operated at 37 mK in a dilution refrig- 
erator in order to minimize the thermal population of the 
excited state. We first characterized each junction sepa- 
rately by performing spectroscopy and Rabi experiments 
(not shown). The junction T\ times are approximately 
3 ns, compatible with the RC time expected for the ~ 
45 pF junctions, where R is an environmental impedance 
of order 100 SI. In Fig. [TJc) we have used the measured 
junction parameters to simulate efficiency r\ for one of 
these junctions as a function of coherent drive frequency 
and detector bias for wr/27t = 200 MHz. The experi- 
mental data shown in Fig. [TJcl) are in good qualitative 
agreement with the model. For a bias point correspond- 
ing to a resonance frequency w 10 /27r = 3.8 GHz, the peak 
detection efficiency is 0.7; the efficiency exceeds 0.6 over a 
band larger than 450 MHz. The appreciable detection ef- 
ficiency and bandwidth of these detectors are well-suited 
to cQED experiments involving high quality on-chip cav- 
ities. 

Next, we examined correlations between the two de- 
tectors for different types of microwave source. The mi- 
crowaves were coupled to the junctions via an on-chip 
power divider with coupling capacitors at each output 
port (Fig. |2ja)). With microwaves applied, we sequen- 
tially pulsed the two junctions into the active state for 
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5 ns, varying the relative delay r of the two measurement 
pulses. The experiments were repeated 10 4 — 10 5 times, 
and we recorded the joint probabilities Poo, Poi, Pio> Pu 
for the two detectors to remain in the supercurrent state 
(outcome "0") or switch to the voltage state (outcome 
"1"). 

In the first experiment, we applied 4 GHz coherent 
photons derived from a commercial microwave generator 
(Fig.[2](a), upper left panel). The microwaves from a com- 
mercial generator constitute a classical electromagnetic 
wave with well-defined phase and amplitude, and provide 
a good approximation to coherent radiation [5J [TH] HH] . 
In the upper panel of Fig. [2jb) we plot the joint switch- 
ing probabilities versus r. We observe no correlations in 
the switching of the two junctions and see that Pn = 
P\Pfi as expected if absorption of a photon is a homo- 
geneous Poisson process. In the second experiment, we 
applied thermal microwaves to the double-junction de- 
tector. The thermal source was realized by amplifying 
and filtering the Johnson noise from a resistor at room 
temperature and upconverting the noise to our 4 GHz 
detection frequency using a microwave mixer (Fig. [2^a) , 
lower left panel). The results of irradiation with the ther- 
mal source are shown in the lower panel of Fig. [2jb) . For 
short relative delay, we observe a clear enhancement of 
Poo and Pu, accompanied by a suppression of Poi and 
Pio- The switching probabilities of the individual junc- 
tions P X A = Pn +Pi and Pf = Pn +P i do not depend 
on the relative delay of the two measurements. 

These experiments constitute a microwave-frequency 
realization of the classic Hanbury Brown and Twiss 
(HBT) experiment |21j , in which temporal correlations in 
the intensity fluctuations at two detectors yield informa- 
tion about the statistical nature of the photon source. We 
define an enhancement factor Ru(t) = Pii(t)/Ph(t — ► 
oo) = Pn(r)/P ] A P 13 . In the standard quantum optics 
treatment of photon counting where the detector is un- 
derstood to interact weakly with the photon field, one 
can show that Rn(r) is equivalent to the second-order 
quantum coherence of the field g^ 2 \r) [7J. For a coher- 
ent state, the photon arrival times are uncorrelated and 
one finds that g^ = 1 for all relative delays, while for a 
thermal state one finds that attains a value of 2 for 
t = and decays to 1 on a timescale set by the coherence 
time of the source. The factor of 2 enhancement of g^ 
is a manifestation of the Bose statistics of the thermal 
source: one observes enhanced probability for measuring 
two (or zero) photons relative to the result obtained for 
a homogeneous Poisson process within a timescale set by 
the characteristic fluctuations of the source. The phe- 
nomenon of photon bunching has been well established 
at optical frequencies by numerous experiments. Clear 
evidence of thermal photon bunching at microwave fre- 
quencies was observed in the cross spectrum of Johnson 
noise measured with two separate chains of linear am- 
plifiers and square-law detectors |19j . Recently, linear 



amplification followed by quadrature detection was em- 
ployed to probe microwave photon antibunching of the 
radiation emitted by a single photon source [22l [23] ; the 
same technique was used to obtain g( 2 '(r) for a con- 
tinuously pumped single photon source and for thermal 
and coherent microwaves |24j . In other work, a protocol 
was developed to perform state reconstruction of propa- 
gating microwave fields from the cross correlations of a 
dual-path amplification setup [25] . and researchers have 
demonstrated state reconstruction of itinerant squeezed 
[2"5] and single photon [57] microwave fields. In most 
of these cases, the sought-after signal is overwhelmed 
by tens of photons of added noise from the microwave 
preamplifier, and extensive signal averaging and postpro- 
cessing of the data are required to extract even lower- 
order temporal correlations, although the authors of [21)] 
employed a Josephson paramp to reduce the added noise 
significantly. By contrast, our Josephson photon count- 
ing scheme provides direct access to temporal correlations 
with a minimal number of preparations of the photon 
state and minimal signal averaging. To our knowledge, 
the data presented here constitute the first observation 
of microwave photon bunching in a coincidence counting 
experiment. 

In Fig. [3] we examine the temporal scale of the pho- 
ton bunching. Here we have varied the bandwidth of 
our thermal source by adjusting the low-pass knee of our 
noise source prior to upconversion. The noise power spec- 
tra are plotted in Fig. [3](a), and the time-domain corre- 
lation data are plotted in Fig. [3|b). We observe that 
the enhancement of Pn decays on a timescale set by the 
inverse of the noise bandwidth, as expected for a multi- 
mode thermal source. 

Finally, we have examined the magnitude of the Pn en- 
hancement as we varied P 1 ' by tuning the input power 
of the thermal source. At low powers, we observe an Pn 
at zero relative delay that is much greater than the factor 
2 expected from the standard quantum optics treatment 
of photon counting. In Fig. [4] we plot the zero-delay en- 
hancement factor Pn(r = 0) versus the single detector 
switching probability \J P^pp . For fixed detector bias 
points, we find that P n scales roug hly as 1/VpApf 
over a broad range of input power, attaining values as 
high as 50 at low power. We have performed experimen- 
tal checks to rule out excess noise from our amplifiers or 
mixers as the source of the enhanced correlated switch- 
ing [28]. While we do not have a definitive explanation 
for the Pn enhancement at this time, there are several 
possibilities. Our strongly coupled Josephson counter in- 
teracts coherently with the input microwave field. In a 
simple model where two detectors are strongly coupled to 
a single thermal mode, one can show that it is possible to 
obtain P n far in excess of 2 for interaction times of order 
half the vacuum Rabi period. We further note that very 
small deviations from a thermal state can dramatically 
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FIG. 3. Temporal scale of photon bunching, (a) Noise power 
spectra of thermal sources with bandwidth 2 MHz (cross sym- 
bol), 4 MHz (dashed line), and 11 MHz (solid line), (b) The 
corresponding joint probabilities Pn as a function of relative 
delay r. The Pn enhancement decays on a timescale set by 
the coherence time of the source, which scales inversely with 
source bandwidth. 



affect the value of Rn. For example, if the occupation 
probability of the n = 2 Fock state is increased by a small 
amount e, then in the limit of small Pi = \J P-^PP, we 
find Rn = er 1 +0(Pi): counterintuitively, the enhance- 
ment of i?n is not bounded by 2 but grows inversely 
with e. We have performed experiments to rule out the 
possibility of an excess of the n — 2 Fock state due to 
downconversion in the detector chip itself. However, we 
cannot rule out the possibility that our noise source gen- 
erates photons with statistics that deviate slightly from a 
thermal state; we plan to investigate this topic in future 
work. 

It should be straightforward to scale this detection 
scheme to tens of parallelized Josephson junctions to en- 
able number- resolved counting of microwave photons. A 
multiplexed Josephson photon counter could be used to 
measure the quantum state of a microwave cavity state 
in the Fock basis. For a measurement time of 10 ns per 
junction and a cavity relaxation time of approximately 
1 /is, tens of junctions can interact with the cavity be- 
fore appreciable decay occurs. Simulations indicate that 
efficient state reconstruction is possible with a handful of 
detectors and a small number of coherent displacements 
of the cavity state. Additionally, a multiplexed Joseph- 
son photon counter would allow measurement of higher- 
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FIG. 4. Pn enhancement Rn = Pu/P^Pf versus 
(P^PP) 1 / 2 for different bias points of the detector; the di- 
rection of the arrow indicates shallower bias points, cor- 
responding to increased dark count rates. For each bias 
point, the probabilities P^' B are tuned by adjusting the in- 
put microwave noise power. The dashed line corresponds to 
Pn = l/^/Pj^PP. At low power, Pn exceeds 2. At lowest 
input powers, switching is dominated by uncorrelated dark 
counts and Pn(r = 0) drops to unity. See text for a detailed 
discussion. 

order correlation functions of the microwave photon field. 
Such a tool could be applied to detailed investigations of 
the statistics of radiation emitted by mesoscopic conduc- 
tors [H] , or to studies of analogue Hawking radiation [55] 
and the dynamical Casimir effect 30 in superconducting 
transmission lines. 
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